A Simplest Swimmer at Low Reynolds Number: Three Linked Spheres 
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We propose a very simple one-dimensional swimmer consisting of three spheres that are linked by 
rigid rods whose lengths can change between two values. With a periodic motion in a non-reciprocal 
fashion, which breaks the time-reversal symmetry as well as the translational symmetry, we show 
that the model device can swim at low Reynolds number. This model system could be used in 
constructing molecular-size machines. 
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The usual swimming mechanism for a human being 
in water involves obtaining a forward momentum from 
the surrounding fluid due to some periodic body motion. 
The fact that the displacement gained in the first half 
period of the cyclic motion is not canceled out by that 
of the second half period is known to be predominantly 
because of the inertial effects Such a mechanism, 
however, does not work in the microscopic world of bi- 
ological objects (such as bacteria), where the effects of 
inertia are not important and the viscous effects domi- 
nate. This case is characterized by very low Reynolds 
number — the dimensionless quantity that measures the 
ratio between the inertial term and the viscous term in 
the hydrodynamical equation of motion 0|. 



Most microscopic biological objects can swim very 
well with velocities of the order of f /im/s, which for 
such micron-sized animals swimming in water, yields 
Reynolds numbers of the order of 10 . In his pioneer- 
ing work, Purcell showed that animals like scallop that 
are equipped with a single hinge cannot swim, using a 
simple opening and closing procedure 0, Q . The reason 
is simple: since the motion is reversible, after finishing 
a cycle the animal will end up being where it initially 
was. He proposed that a non-reciprocal motion, which 
breaks the time-reversal symmetry, is needed to produce 
a net displacement. This will help the animal to pro- 
pel itself, during each cycle, along some direction that 
is preferred by the symmetry of the system and the mo- 
tion [31 . Despite the simplicity of the geometry that he 
suggested — three rigid rods connected with two hinges — 
quantitative analysis of the "Purcell swimmer" had not 
been performed until very recently, due to the complex- 
ity of the Stokes equation in the specific geometry 
We note that there are other related swimming mecha- 
nisms at the mesoscopic scale in the context of molecular 
machines Q- 

Here we use Purcell's original idea and introduce a very 
simple and experimentally accessible model system that 
can swim using proposed periodic internal motions. The 
swimmer consists of three hard spheres that are linked 
through two arms, and has the advantage that the de- 
tails of the hydrodynamic interactions, as well as the 
swimming velocity and direction, can be worked out with 



FIG. 1: Three linked spheres connected by two arms of neg- 
ligible thickness. 



great ease, as compared to the case of the Purcell swim- 
mer. 

The model swimmer that we are proposing is shown 
in Fig. ^ an d consists of three spheres with radius R 
that are connected by rigid slender arms aligned along 
the x direction. The spheres are assumed to be floating 
in a highly viscous fluid with viscosity /z. There are two 
internal engines on the middle sphere (sphere number 1), 
which act as internal active elements responsible for mak- 
ing a non-reciprocal motion that is needed to propel the 
whole system. We consider the initial state of the sys- 
tem such that the spheres number 2 and 3 are in equal 
distance D from the middle sphere. We divide a com- 
plete cycle of the non-reciprocal motion into four parts 
as below (see Fig. 

(a) In the first step of the motion, the right arm has fixed 
length, and the length of the left arm is decreased with 
a constant relative velocity W, using one of the internal 
engines in the middle sphere. We denote the relative 
displacement of the spheres 1 and 2 in this stage by e. 

(b) As the second step, the left arm is fixed and the right 
arm decreases its length with the same constant relative 
velocity W as before. The relative displacement of the 
spheres 1 and 3 is again e, like the previous stage. 

(c) During this step, while the right arm is kept fixed, 
the left arm increase its length with the same relative 
velocity W to reach its original length D. 

(d) Finally, in the last step the left arm is kept fixed 
and the right arm elongates to its original length with 
the same constant velocity W. The system is now in its 
original internal configuration. 

As can be seen from Fig. [21 the above 4-stage cycle 
is not invariant under time reversal, and we can thus 
expect a net translation upon completing a full cycle. To 
obtain a net translational motion, the above cycle can be 
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FIG. 3: An auxiliary (fictitious) movement in which the right 
arm has a constant length S while the left arm changes its 
length from D to D — e. During this movement the middle 
sphere will be displaced by an amount Af (8). 



FIG. 2: Complete cycle of the proposed non-reciprocal mo- 
tion of the swimmer, which is composed of 4 consecutive time- 
reversal breaking stages: (a) the left arm decreases its length 
with the constant relative velocity W, (b) the right arm de- 
creases its length with the same velocity, (c) the left arm opens 
up to its original length, and finally, (d) the right arm elon- 
gates to its original size. By completing the cycle the whole 
system is displaced to the right side by an amount A. 



repeated continuously. 

The general equation that describe the hydrodynam- 
ics of low Reynolds number flow is the Stokes equation 
for the velocity field u, subject to the incompressibility 
condition: 



/xV 2 u 
V • u = 



o, 



0, 



(1) 

(2) 



where p represents the pressure field in the medium. 

Assuming that the spheres are moving inside the fluid 
with velocity vectors , with the index i denoting the la- 
bels of the spheres, the description of the system involves 
solving the above-mentioned equations with zero veloc- 
ity boundary condition at infinity and no-slip boundary 
conditions on the spheres, which implies 



U|r on the i'th sphere 



V,; 



The variables that are necessary to determine the dynam- 
ics of the spheres are their velocities and the forces 
Fi acting on them. By solving the above governing equa- 
tions, we will be able to obtain the fluid velocity in the 
medium, and hence, the corresponding stress tensor that 
will give us the required forces on the spheres. It is a 
simple observation that because of the linearity of the 
governing equations, and the linearity of the stress ten- 
sor with respect to the velocity field, one can generally 
expect a relation of the form: 



(4) 



where the symmetric Oseen tensor Hij depends on vis- 
cosity, the geometry of the bodies immersed in it (in our 



case the spheres), and their relative orientations. Includ- 
ing the condition that there are no external forces such 
as gravity, the system of spheres should be force-free: 



i=i 



0. 



(5) 



Since we are only interested in the dynamics of the 
spheres, we can equivalently solve the set of Eqs. and 
(J5J, instead of Eqs. QJ and @. To solve these equations 
we need to know the form of the Oseen's tensor. Let us 
consider a coordinate system in which the position vector 
of the i'th sphere is x, and the separation between the 
i'th and the j'th spheres will be Xy- = Xj — Xj, with a unit 
vector n in this direction. General symmetry considera- 
tions will allow us to write the hydrodynamic interaction 
tensor in the following form [3, |9| : 
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[ J 4 ii (A)nn + S ij (A)(I-nn) 



(6) 



where we use the dimensionless quantity A = R/xij. As- 
suming that the separations between the spheres are suf- 
ficiently larger than their sizes, we can write a perturba- 
tion expansion for the symmetric coefficients Aij and By 



(3) in powers of A, which reads 



and 



A, 



1 + 0(A 4 ), 

§Ay + 0(A 3 ), 



1 + 0(A 4 ), 
f X i3 + 0(A 3 ), 



i = J 



» = J 



(7) 



(8) 



to the leading order. Note that for simplicity we are con- 
sidering only the translational motion for the spheres, al- 
though the effects of rotational motion can be taken into 
account in a similar way. Moreover, extra simplification 
comes from the one-dimensional nature of the motion, 
such that the tensorial structure of the hydrodynamic 
interactions plays no important role in the dynamics. 

To analyze the motion of the system during one com- 
plete period of the non-reciprocal cycle, we introduce an 
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auxiliary move as shown in Fig. [3J During this mo- 
tion, the right arm has a constant length 5 while the left 
arm changes its length from D to D — e with the con- 
stant velocity W. Using symmetry arguments, we can 
relate all the four steps in the non-reciprocal cycle to the 
above move, as follows: The step (a) corresponds to the 
auxiliary move by setting S = D. By applying a reflec- 
tion transformation with respect to any point on the x 
axis, the step (b) can be mapped onto the auxiliary move 
with 5 = D — e. The step (c) is obtained by applying a 
time-reversal transformation on the auxiliary move, with 
5 = D — e. Finally, the step (d) is obtained by applying 
a reflection transformation [as in step (b)] followed by a 
time-reversal transformation, with S = D. 

To obtain the net displacement of the middle sphere 
in the real problem, it is thus enough to solve the dy- 
namical equation for a single auxiliary movement. If we 
define the net displacement of the middle sphere dur- 
ing the auxiliary step by Af(<5), then by considering the 
above arguments we can calculate the total displacement 
A of the real system through a complete cycle as 

A = 2[A f (D) - At (D - e)]. (9) 

Noting that the above displacement takes place during 
the time of a complete cycle that is 4e/W, then we can 
calculate the average swimming velocity for the swimmer 
as 

v s = W ^. (10) 

As an example, we have numerically calculated the dis- 
placements for the case of D = 10 R and e = 4 R. During 
the first stroke the middle sphere swims in the — x direc- 
tion by an amount 1.35 R, while in the second stroke 
it swims a distance of 1.44 R in the positive direction. 
The third stroke then causes a continuation of the mo- 
tion in the positive direction for a distance of 1.44 R, and 
finally, the fourth stroke takes it back by a distance of 
1.35 R. At the end of this cycle the sphere is displaced 
by a net amount of 0.16 R in the positive direction, as 
shown in Fig. [21 Figure 01 shows the total displacement 
of the middle sphere during a complete non-reciprocal 
cycle as a function of the internal relative displacement 
e for D = 10 R. 

In the limit of small internal deformation of the swim- 
mer, we can expand all the quantities in terms of e/D, 
and calculate the swimming velocity in a perturbative 
series. To the leading order, we find 

v.- 17* (£)(£)■. <n> 

The above result shows that the scale of the swimming 
velocity is set by the typical velocity of the internal mo- 
tion. Moreover, the swimming appears to be a quadratic 
effect with respect to the small deformations in the sys- 
tem. These two characteristics are general, as can be 
seen in other swimmers at low Reynolds number. 
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FIG. 4: Dimensionless displacement of the swimmer in a com- 
plete cycle as a function of the dimensionless relative displace- 
ment between neighboring spheres. 



In our simple prescription for the non-reciprocal mo- 
tion we have assumed that in each step one sphere moves 
with respect to the middle one and the other one is kept 
in relative constant distance. This assumption has been 
made for simplicity, and we can imagine a more gen- 
eral continuous motion of the spheres with respect to the 
middle one. The only necessary condition is to break 
the time-reversal symmetry and obtain a non-reciprocal 
motion. The general requirement in the continuous mo- 
tion will be the existence of a nonzero phase difference 
between the continuous periodic motions of the left and 
right spheres with respect to the middle one. The swim- 
ming velocity in this case will also have the two general 
characteristics mentioned above. 

The generalized case of many spheres that are coupled 
to each others on a regular one- or two-dimensional lat- 
tice is a suitable microscopic model for an extensible flat 
filament or membrane. In these cases, the relative "in- 
plane" motion of the neighboring spheres with respect 
to each other will cause the whole system to swim. The 
direction of motion and the swimming velocity can be 
predicted using the simple proposed description. The in- 
ternal relative motion of the system that can cause swim- 
ming is of the general form of a traveling wave on the 
position of the spheres in some direction. Traveling wave 
is the simplest non-reciprocal motion which breaks the 
time-reversal and translational symmetries. The swim- 
ming velocity will be opposite to the direction of the 
wave-vector, and proportional to the phase velocity of 
the traveling wave. Such kind of motion has been con- 
sidered for a spherical membrane, and shown to have the 
same characteristics as suggested above [To| . 

In conclusion we have introduced a very simple swim- 
mer and calculated its swimming velocity. The swimmer 
uses some periodic internal motion to propel itself under 
low Reynolds number conditions. The advantage of this 
model, as compared to previously known model swim- 
mers, is that the analysis of the hydrodynamics problem 
can be performed considerably more easily. The model 
swimmer could be used in making molecular-size ma- 
chines with controllable motion. 

We have benefited throughout this work from fruitful 
discussions with A. Ajdari, F. Julicher, K. Kruse, and H. 
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